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The dispersive interaction of atoms and a far-detuned light field allows nondestructive imaging of 
the density oscillations in Bose-Einstein condensates. Starting from a ground state condensate, we 
investigate how the measurement back action leads to squeezing and entanglement of the quantized 
density oscillations. In particular, we show that properly timed, stroboscopic imaging and feedback 
can be used to selectively address specific eigenmodes and avoid excitation of non-targeted modes 
of the system. 
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Nonclassical, squeezed and entangled states of light 
provided by nonlinear optical components are important 
ingredients in quantum metrology and communication 
HE]. Production of analogous states of matter have been 
a long quest in physics with implementations demon¬ 
strated in atomic ensembles, as well as in superconduct¬ 
ing and nanomechanical devices. Relying on nonlineari¬ 
ties due to atomic interactions, ultracold atoms have been 
prepared in two-mode entangled states of split atomic 
clouds sill and internal state components [bUTO]. Fur¬ 
thermore, it has been possible to prepare states that wit¬ 
ness fundamental quantum properties through violation 
of inequalities obeyed by classical models [nHE]. 

Using dispersive light-matter interactions and mea¬ 
surement back action, room temperature vapour exper¬ 
iments have demonstrated squeezed HZ] and entangled 
m states, quantum teleportation m, and a quantum 
memory for light m- Similar experiments with interact¬ 
ing cold atoms and light fields have shown great progress 
[ 2 TII 24 ] . and numerous proposals exist [251127] to manip¬ 
ulate the collective quantum states of ultracold atoms, 
while the full exploitation of the quantum mechanical 
nature of the interaction and measurements is yet to be 
realised. 

In this Letter we develop a theoretical treatment of 
Bose-Einstein condensate (BEG) dynamics due to the 
dispersive interaction with an optical probe field, sub¬ 
sequently measured, and feedback [Fig. [^a)]. Specifi¬ 
cally, we investigate the selective preparation of squeezed 
and entangled states of the multimodal excitations of the 
BEG. Due to the temporal evolution, the local atomic 
density is not a quantum non-demolition (QND) vari¬ 
able [28U3Q] , but by applying stroboscopic measurements 
at selected times, we can address effective QND variables 
of selected eigenmodes of the BEG dynamics. 

We consider a ID ultracold Bose gas m harmonically 
confined with axial (radial) trapping frequency uJx{D) 
the length scale where radially tight 

confinement restricts the low energy excitations to ax¬ 
ial motion. The BEG ground state wave function /q {x) 
(taken to be real) is given by the ID Gross-Pitaevskii 


equation 

[Hi + gi-Dni,{x)]fQ (x) = h/q (x), 

where the EEC meanfield density is no(x) = V/g (a;), and 
the chemical potential, /i, enforces the BEG population 
to the total number of atoms, N. Hi is the single atom 
Hamiltonian and the ID interaction strength is gm = 
with the 5- wave scattering length 
The elementary excitations of the BEG are the collec¬ 
tive center-of-mass, breathing, and higher order modes, 
where each mode, jf, constitutes a quantum harmonic os¬ 
cillator degree of freedom with (dimensionless) quadra¬ 
ture observables Xj{t) and Pj(t). Their frequencies, ujj, 
and wavefunctions, ff{x), are found by solution of the 
coupled Bogoliubov-de Gennes equations. 
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where C± = i7i — /i + (2 ± 1 )^id^o (^) • These eigenmodes 
provide an expansion of the probed atomic density, 

h{x,t) = n(i{x) + 2y^^iJno{x)f~{x)xj{t) + 0{N°), (1) 
j 

embodying the dynamic quantum fluctuations about the 
BEG meanfield. As we shall see, stroboscopically probing 
the density enables mode selective squeezing [red line in 
Fig. [^b)] and entanglement [heralded by nonzero off- 
diagonal elements in Fig. [^c,ii)], and reinitialization of 
the modes to oscillator ground states [end of Fig. [^b)]. 

The dispersive light-matter interaction is characterized 
by the coupling constant k = with the atomic 

depumping rate g, and optical depth on-resonance, do 
[32l [ 33 ] . As we image the spatial density 0 by optical 
phase shift measurements [Fig. [^a)], the light field de¬ 
tection in a certain pixel is sensitive to (a combination 
of) the atomic variables Xj{t). This is described by the 
atomic-pixel mode overlap integrals, 

fjd = jdx jdx'Kiix - x')fQ {x')f~ (x'), 
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FIG. 1. (Colour online) (a) A BEG imprints a phase shift on a planar coherent light field depending on the atomic density along 
the BEG axis. The phase shifts are spatially detected by an encapsulating (/l > lx) array of homodyne detectors (pixels) of 
widths Id = /x/10. (b) The atomic modes can be selectively addressed and squeezed, illustrated here, where first the 3^^ mode, 
and subsequently, modes 1 and 5 are squeezed, (c) The absolute value of the covariance matrix elements are shown for the first 
three odd modes. Continuous probing (i) squeezes and correlates a swathe of atomic modes, while stroboscopic probing (ii) can 
generate correlations between only selected modes. The same level of entanglement of modes 1 and 3 is achieved in (i) and (ii), 
while the strong squeezing and the addressing of 5^^ in (i) is absent in (ii). The simulations are performed for probing the D 2 
line (cr+ polarised light) of 1000 ®^Rb atoms in the \F,Fz) = |2, 2) state with Ux — ^tt x 150Hz, u± = lOOct;^, and = 2fiujx- 


where dx denotes integration over the domain of the 
pixel. The convolution with the measurement ker- 

nel [34l [35], JCa{x) = Jdke e*^^/27r, accounts 

for the resolution limit along the BEC axis associated 
with the diffraction of light propagating over distances 
^ l± through the BEC. By smearing out spatial features 
smaller than /r = (/xA) 2 , where A is the light wave¬ 
length, it prevents the addressing of the higher modes 
with shorter spatial variations (the wavelength of mode 
j scales as ^ Ixjj)- We will also require, for our analy¬ 
sis of the correlations between modes, the atomic mode 
overlap integrals 

fjk = JdxJdxX^{x-x')fo{x)f-{x)fo{x')f-{x'). 

Since the dominant contributions of the light-matter 
interaction are second-order in light and atomic quadra¬ 
tures, and the light field is subject to quadrature mea¬ 
surements, the atomic quantum state can be described by 
a gaussian Wigner function. At all times, the state is fully 
characterised by [R]j = {uj) and [A]jk = cov(z>j, i>/c), 
the first and second moments of the atomic mode vari¬ 
ables, z> = [xi{t),pi{t),X 2 {t),p 2 {t),.. .]^^ respectively. 
Following the methodologies of [36|, the quantum back- 
action of the light field measurements of the p quadrature 
in each pixel [represented by expectation values and ran¬ 
dom Wiener increments, dWd{t)] results in the stochastic 
evolution of the first moments [33] 

dR = — DRdt + AmdW. (2) 

The harmonic rotations in all {xj{t)^pj{t)} phase spaces 
is represented by the block diagonal matrix, D, with 2x2 


blocks [Dio- = 


0 


-UJj 

0 


Each pixel probes a linear 


combination of modes, represented by a rectangular ma- 

0 2Kfjd 


0 0 


, and owing 


trix of 2X2 blocks [m ]— ~y ^ 

to the correlations among the modes (A), the measure¬ 
ment results, dW = [0, dWi{t),0, dW 2 {t), .. .]^, affect the 
modes in a correlated manner. 

The covariance matrix (second moments) evolve as [33| 


A = E - DA - AD^ - AMA, 

where M = mm^ and the matrix E with 2 

0 0 


[^]jk — 


0 


jk 


(3) 

X 2 blocks, 
A evolves in a manner indepen¬ 


dent of the measurement outcomes. While this implies 
we can deterministically assess the squeezing and entan¬ 
glement properties of the system, we recall the random 
measurement results determine the mean displacements, 
about which, the squeezing and entanglement occurs. We 
shall now first consider the covariance matrix solutions 
of Eq. focusing on the squeezing and entanglement 
generation, and later return to the stochastic evolution 
of the mean displacements described by Eq. <§■ 

Probing with a constant optical field strength, 
squeezes Xj{t) at a rate Uj = while the conju¬ 

gate quadrature Pj{t) is anti-squeezed. Since Xj{t) and 
Pj{t) are coupled by rotation at rate cjj, the squeez¬ 
ing may be effectively suppressed. This is reflected by 
Eq. @ yielding a steady state solution where the Xj{t) 
variance var[xj(t)]ss — [(1 + 4z/|)^ — 1]^/2^/2i?j < 1/2, 
and Uj = Only if appreciable squeezing occurs 

during a fraction of a full phase space rotation, will there 
be squeezing in the long time limit. 
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For a given density, the rate of squeezing can only 
be increase by increasing the atomic depumping rate, 
T]. Alternatively, a squeezed state of an oscillator with 
quadrature Xj{t) = cosujjt sinujjt, where the x^ 
quadrature is squeezed, can still be achieved by apply¬ 
ing a temporally modulated field m- Here, we imple¬ 
ment stroboscopic probing with a sequence of constant 
intensity pulses determined by one or more frequencies, 
zji. The timings, t/, and durations, r/, of the pulses 
are determined by ensuring, for each frequency, zjit is 
within lS.pl2 from a multiple of 27r. Probing at the sin¬ 
gle frequency, w = 2co’j, amounts to a train of n pulses 
centred on times U = [0, Tr/cjj, 27r/co’j,..., nTr/cjj] with 
identical durations, ri = r = Apjw. Such a field ad¬ 
dresses, and squeezes, the quadrature, while avoid¬ 
ing its anti-squeezing at intermediate times. Thus, en¬ 
abling squeezing well beyond the continuous case for 
the same strength, k.. For small Ap, the effective 
QND probing of x^ is expected, and hence, the variance 
var[x^]QND = 1/2{1 P 2i/jTT) [38] for the accumulated 
probing time tt = nr. This is observed (red line) in 
Fig-Sb), where the first quarter of the trace squeezes the 
3^^ mode = 100co’a,/27r, zu = 2uj^ and A(^/27r = 0.05). 

The atomic multimodal system is far from the sin¬ 
gle mode squeezing picture. However, if no other mode, 
nor coupled correlation, has a rational frequency ratio to 
cjj, the pulse train only addresses, and squeezes, the 
mode. This modal selectivity is illustrated in Fig. [^b), 
where the variances of modes 1 and 5 are essentially 
unaffected during the probing of the 3^*^ mode. Here, 
the atomic interactions provide an irregular spectrum of 
the lowest frequency modes m, allowing separate ad¬ 
dressing without crosstalk. In light of this, evaluating 
the performance of squeezing and entangling operations, 
we will later assess the selectivity within the multimodal 
j = 1,3,5 subspace, as well as, this subspace’s isolation 
from the rest of the system. 

The simultaneous squeezing of two modes, j and k, 
is achieved with the two frequencies, wi = 2ujj and 
W 2 = 2co’/c, as featured in Fig. [^b) for modes 1 and 5 
after ujxt = 257r {Ap/2r: = 0.1). The probing sequence is 
now out-of-phase with the 3^*^ mode, and hence, its prior 
squeezing is progressively lost returning to initial vacuum 
values. The final 5% of the trace demonstrates reinitial¬ 
ization of all modes to vacuum fluctuations by switch¬ 
ing to weak continuous wave probing = 50ci;a^/27r). 
The preservation of vacuum fluctuations of nontargeted 
modes, loss of squeezing of the 3^^ mode, and reinitial¬ 
ization of all modes to vacuum fluctuations, featured in 
Fig. Bb), are demonstrations of the interplay between 
atomic dynamics, measurement strength, and strobo¬ 
scopic probing, enacting a quantum eraser for a subset 
of, or all, atomic modes. 

To further assess the performance of the operations, we 
study squeezing of the 3^^ mode in Fig. (similar results 
are found for other modes and for the joint squeezing 



FIG. 2. (Color online) The squeezing of the 3’^'^ mode after 
stroboscopically probing (/^^ = 100ct;x/27r) for 100 trap peri¬ 
ods is shown, while the temporal evolution of the quadrature 
variance, var[:r 3 ], is illustrated in the inset. For small A(j) 
(pulse duration), var[x 3 ] follows the QND result (red line, see 
text), while for larger A0, the squeezing is faster but subop- 
timal. Excellent selectivity, Hellinger distance, Dq-c ~ 0, and 
little crosstalk, Purity, P ^ 1, are observed for small A(j). The 
results for the noninteracting atomic system (dotted lines) are 
shown for comparison. 

or entangling of pairs of modes). The rate of squeez¬ 
ing is determined by A0, as the accumulated probing 
time Tt oc Acj) (inset). However, only for smaller A0, 
the x-quadrature is probed in a QND fashion (red line), 
while the squeezing for larger A0 is suboptimal as a re¬ 
sult of the inadvertent probing of the p-component. The 
figure also addresses crosstalk between the j = 1,3,5 
subspace and its relative complement through the pu¬ 
rity, P, of the reduced system density matrix, p. The 
Hellinger distance [39], = Tr[p 2 — (pc^)2]^/2, quan¬ 

tifies the selectivity within the subset of modes. The 
desired state pd assumes identical R to p, but with all 
blocks [A]j/c, except j = k = 3^ replaced by their ini¬ 
tial vacuum values. Excellent selectivity {D'^ 0) and 

little crosstalk (P 1) are observed for small A0. For 
comparison (dotted lines), a noninteracting atomic sys¬ 
tem demonstrates similar squeezing, however, the corre¬ 
sponding linear spectrum ujj = ujxj results in significant 
crosstalk and poor mode selectivity. 

A probe transmitted through two atomic media reveals 
information about their collective rather than individual 
properties and may hence lead to their mutual entan¬ 
glement [38]. Similarly, we may entangle two modes, j 
and /c, of a single BEG, by probing the density oscillation 
amplitude in a manner that does not discriminate contri¬ 
butions from the individual modes. The modes’ respec¬ 
tive spatial signatures, ^ /o ix)fj (x) and ~ (a;)/^. (x), 

must be indistinct {fjf. ^ 0). Partial temporal distin- 
guishability, owing to different oscillation frequencies, ujj 
and cj/c, is avoided by stroboscopically probing with a 
train of pulses at the single frequency, w = cjj + cj/c. 
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FIG. 3. (Color online) (a) A pulse sequence of 100 trap peri¬ 
ods with zui = cji + 6 ^ 3 , /\lpI2'k — 0.03, and h? — 30ct;cc/27r, 
entangles modes 1 and 3. The logarithmic negativity E 12 , 
slightly out-performs the corresponding QND result, F^qnd- 
The selectivity within the subspace is excellent {Du ~ 0), 
while the subspace’s purity P ~ 97%. (b) The mode func¬ 
tions (not to scale) of the condensate, fo{x), entangled modes, 
/i“ 3 (x), and the overlap, ~ fo{x)fr {x)U (x)- 



FIG. 4. (Color online) A stochastically evolving trajectory of 
{x 3 {t)) and {p 3 {t)), and the standard deviations of 1000 tra¬ 
jectories, cr(^)(t) and cr(p 3 )(t), corresponding to the simula¬ 
tion in Fig. 0 with A0/27r = 0.15. The grey regions indicate 
when the measurement occurs. Without feedback (a), the 
standard deviations of the traject ories can be approximated 
by cr^^ 3 ^(t) = ^{p 3 ){t) — \/r' 3 TT /2 (red line). With feedback 
(b), we demonstrate the desired damping of mean values. 


Analogous to the case of squeezing, it also allows selec¬ 
tive addressing of modes as illustrated in Fig. Be). The 
absolute value of the covariance matrix elements of the 
j = 1,3,5 subspace is shown for (i) the steady state of 
continuous probing {k‘^ = 1000co’a;/27r) and (ii) strobo¬ 
scopic probing for 100 trap periods (/^^ = 4co’a,/27r and 
Aipl27r = 0.03). Although both cases reach the same 
level of entanglement of modes 1 and 3, the continuous 
probing addressing a swathe of modes, causes an addi¬ 
tional squeezing, and requires a larger strength, k. 

Bipartite entanglement can be quantified by the log¬ 
arithmic negativity Ej^ = log 2 ||p^^||xr 02]’ where the 
partial transpose and trace norm are denoted, Tp, and 
ll'llrp^, respectively. Figure |^a) shows the temporal build 
up of entanglement between modes 1 and 3, and the com¬ 
parison to an effective QND probing that converges for 
^ 00 to the asymptotic limit, ^qnd ^ log 4 [i^]- 
The spatial distinguishability of the mode functions, fea¬ 
tured in Fig. ib), remains the limiting factor, and is 
parameterised hy /3 = \fjk\/[fjjfkk\^’ Deviations from 
^QND attributed to coupling to other modes with 
commensurate correlation frequencies [here, uJi P 
2{uJi + cj?)] through as signalled by a small loss of 
purity (P ^ 97%). The selectivity within the subspace 
is excellent {Du — 0). 

We have so far demonstrated that stroboscopic probing 
permits squeezing and entanglement of particular reso¬ 
nant modes. As illustrated in Fig.j^a), due to the prob¬ 
ing backaction, the modes are also subject to random dis¬ 
placements, governed by the diffusion terms in Eq. <§■ 
The Gaussian multi-mode quantum state is, indeed, sub¬ 
ject to such coherent displacement of all modes. This cor¬ 
responds to a modification of the Gross-Pitaevskii mean 


field wave function, and it is suggestive that modula¬ 
tion of the trapping potential can be used as a feedback 
mechanism to recover and maintain the original shape 
fo{x). A micro-mirror array [41] or spatial light modu¬ 
lator [42] can provide such an adaptive lightshift poten¬ 
tial for the atoms. Choosing the potential strength with 
single atom Hamiltonian, Hf = hk{x){pk{t)), where 


hk{x) = {x)5jk/ y/no(x), we can selectively address 

the displacement of the mode. The corresponding 

action on the mode is the addition of ^ to 

[ 0 2ujj 

the 2x2 block [D]jj in Eq. ([^. Due to the degenerate 
eigenvalue —ujj , the deterministic part of the evolution of 
Xj {t) and pj {t) becomes critically damped and suppresses 
the displacements caused by the dW-terms in Eq. 

We have simulated the full multi-mode dynamics, and 
Eig. [^b) demonstrates the displacement of the squeezed 
mode in Eig. is successfully suppressed. 

The feedback demonstrated in Eig. [^ is only performed 
for a single mode (although not limited to). In the ef¬ 
fective QND examples of single mode squeezing (e.g., 
Acj) = 0.01 in Eig.[^ and bipartite entanglement (Eig.[^, 
the coherent excitation of mean values corresponds to an 
^ 8% population A^nc outside of the BEC mode. If neces¬ 
sary, more elaborate can be investigated to minimise 
the excitation of multiple modes, and N^c iiiay be re¬ 
duced by different trapping geometries, where the mea¬ 
surement kernel 1C a {x) may yield less coherent excitation 
of the irrelevant modes. 


In conclusion, we have demonstrated the quantum con¬ 
trol of a matter-wave system via spatially resolved opti¬ 
cal probing. Using stroboscopic probing and feedback, 
we can address and correlate effective QND observables 
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of selected density modes of a BEC, while preserving 
the initial vacuum fluctuations of nontargeted modes. 
The performance of our proposal is ultimately limited by 
residual atomic depumping [43] and by the depletion of 
the BEC ground state when many modes become signif¬ 
icantly squeezed. With current experimental parameters 
[44h46] . however, a small alkali BEC should undergo suffi¬ 
cient squeezing and entanglement to allow demonstration 
of a few-mode memory for scalable quantum information 
processing. 
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